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We test the predictive power of first-oder reversal curve (FORC) diagrams using simulations of 
random magnets. In particular, we compute a histogram of the switching fields of the underlying 
microscopic switching units along the major hysteresis loop, and compare to the corresponding 
FORC diagram. We find qualitative agreement between the switching-field histogram and the 
FORC diagram, yet differences are noticeable. We discuss possible sources for these differences and 
present results for frustrated systems where the discrepancies are more pronounced. 
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INTRODUCTION 

The conventional methods to characterize mag- 

netic interactions in hysteretic systems, such as the 8M 
method 0,0]) utilize isothermal remanent magnetization 
(IRM) and dc demagnetization remanence (DCD) curves 
based on the Wohlfarth relation [(| . Recently, FORC di- 
agrams 0, @ have been introduced to study hysteretic 
systems. Their extreme sensitivity has helped to "finger- 
print" several experimental systems as well as theoretical 
models ranging from geological samples and recording 
media to paradigmatic models of random magnets and 
spin glasses §• 

In this work we perform numerical simulations of ran- 
dom magnets (and spin glasses) in order to test the pre- 
dictive power of FORC diagrams by comparing to a his- 
togram of up- and down-switching fields of the underlying 
switching units along the major hysteresis loop. 

The aforementioned re-parametrization of the major 
hysteresis loop (switching-field histogram) displays the 
information carried by the major loop in a more com- 
prehensive way and provides a good comparison to the 
FORC diagram. We find, that the major-loop behav- 
ior predicts the minor-loop behavior captured by FORC 
diagrams well. We present a comparison of both distri- 
butions and discuss some differences between them. We 
argue that switching-field histograms arc useful to study 
hysteretic systems in more detail than with conventional 
methods due to their simplicity and ease to compute. 

MODEL & ALGORITHM 

The Hamiltonian of the random-field Ising model 
(RFIM) is given by 

h = J » s i s i - Yl hiSi -Hj2 Si - w 

Here Si = ±1 are Ising spins on a square lattice of size 
N = L 3 in three dimensions with periodic boundary con- 
ditions. The interactions between the spins are uniform 



1) and nearest-neighbor, and H represents the 



externally applied field. Disorder is introduced into the 
model by coupling the spins to site-dependent random 
fields hi drawn from a Gaussian distribution with zero 
mean and standard deviation or. 

We simulate the zero-temperature dynamics of the 
RFIM by changing the external field H in small steps 
starting from positive saturation. After each field step 
we compute the local field fi of each spin: 



fi — Y^ JvSj — a — hi 



(2) 



A spin is unstable if it points opposite to its local field, 
i.e., fi ■ Si < 0. We then flip a randomly chosen unstable 
spin and update the local fields at neighboring sites. This 
procedure is repeated until all spins are stable. 

For the rest of this work we set L = 50 (N = 125000 
spins) and or = 5.0, unless otherwise specified. The 
different figures are calculated by averaging over 5000 
disorder realizations in order to reduce finite-size effects. 



FORC DIAGRAMS 

In oder to calculate an FORC diagram, a family of 
First Order Reversal Curves (FORCs) with different re- 
versal fields -£Zr is measured, with M(H, £zr) denoting 
the resulting magnetization as a function of the applied 
and reversal fields. Computing the mixed second order 
derivative plLuj 



p(H,H R ) = ~[d 2 M/dHdH R ] 



(3) 



and changing variables to H c = (H — H-g)/2 and = 
(H + Hr)/2, the local coercivity and bias respectively, 
yields the "FORC distribution" p(H h ,H c ). FORC dia- 
grams resemble the commonly known Preisach diagrams 
[llL Il2| , yet they are model-independent and therefore 
more general. 

Figure Q shows an FORC diagram of the RFIM at 
high disorder strength (or = 5.0 > a c nt ~ 2.16) 
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FIG. 1: FORC diagram of the RFIM for disorder a R = 5.0, 
well above the critical disorder 13] in three dimensions. Note 
the pronounced vertical feature at H c = 1 with a wake extend- 
ing to H c = which corresponds to multi-domain nucleation 
in the sample. The dots along the H^-axis are numerical noise 
(no data smoothing). 

Note the vertical ridge at H c « 1 reminiscent of domain- 
wall nucleation 0] • A vertical cross-section of the ridge 
(H c = 1) mirrors the distribution of the applied ran- 
dom fields, because these can be viewed as a distribution 
of random biases acting on the spins when ctr, 3> cr C rit- 
We have tested this in detail by selecting the random 
fields from a box distribution. The resulting FORC di- 
agram is qualitatively similar to the Gaussian case, yet 
a vertical cross-section of the ridge is box-shaped. This 
is not evident by studying the major hysteresis loop for 
different disorder-distribution shapes and illustrates the 
advantages of the FORC method over conventional ap- 
proaches for studying hysteretic systems. 



SWITCHING-FIELD HISTOGRAMS 

In order to test the predictive power of FORC dia- 
grams, we simulate the RFIM with the zero-temperature 
dynamics described in Sec. D and store the up- and down- 
switching fields of the spins along the major hysteresis 
loop. We then create a histogram of the number of 
flipped spins for a given pair of up- (Hi) an( i down- 
switching fields (Hi). By changing the variables to 
the coercivity [H c = (Hi — Hi)/ 2] and bias [H^ = 
(Hi +Hi)/2] of each spin, we obtain a distribution of the 
coercivities and biases of the spins in the system along 
the major hysteresis loop. 

Figurc[21shows the switching-field histogram (SFH) for 
the RFIM. One can see a close resemblance with the cor- 
responding FORC diagram presented in Fig. ^ In order 
to better compare FORC diagram and SFH, in Fig. [3] we 



FIG. 2: Switching-field histogram of the three-dimensional 
RFIM with or = 5.0. Note the close resemblance to the 
FORC diagram presented in Fig.0 Because no derivatives of 
the data are required, the contours are much smoother than 
in the case of an FORC diagram. 
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FIG. 3: Absolute difference between the FORC diagram pre- 
sented in Fig. and the corresponding SFH in Fig. [5] for the 
three-dimensional RFIM with gr = 5.0. The details are dis- 
cussed in the text. 

present the absolute difference between both diagrams. 

Even though the SFH and the FORC diagram of the 
RFIM differ slightly (see Fig. the main character- 
istics representing the underlying physical properties of 
the model are the same (vertical ridge representing multi- 
domain nucleation). It is interesting that a zeroth-order 
reversal curve (the major hysteresis loop) contains possi- 
bly all the necessary information to reconstruct the first 
order reversal curves of the system (the FORC diagram) . 

The differences found between the FORC diagram and 
the SFH could be due to numerical error in the deriva- 
tives of the FORCs because noise is amplified in numer- 
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ical derivatives considerably. In addition, the differences 
could be attributed to either hysteron correlations or the 
failure of the (simple) Preisach picture of hysterons. The 
latter would imply that a generalization of "classical" 
hysterons is required. 

Figure |3 also illustrates how the re-parametrization of 
the major hysteresis loop in terms of an SFH shows more 
details about the microscopic structure of the system. 
The gained information is similar to the information pro- 
vided by an FORC diagram, yet the computation of an 
SFH is considerably faster than calculating an FORC di- 
agrams (generally ~ 10 2 times faster) and involves no 
numerical derivatives of the data, thus reducing numeri- 
cal error. 
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FRUSTRATED SYSTEMS 

As the random-field Ising model is a random magnet 
with no frustration, we also calculate the FORC diagram 
and SFH for the 3D Edwards- Anderson Ising spin glass 
15] (EASG). Due to frustration, a spin can flip more than 
twice along the full hysteresis loop. With the current 
definition of the SFH this is not taken into account and 
differences to an FORC diagram are expected. 

The Hamiltonian of the Edwards- Anderson Ising spin 
glass is given by Eq. JIJ where the Jij are nearest- 
neighbor interactions chosen according to a Gaussian dis- 
tribution with zero mean and standard deviation unity, 
and hi — Vi. H represents the externally applied mag- 
netic field and periodic boundary conditions are applied. 
For the simulations we use the zero-temperature algo- 
rithm presented in Sec. B Frustration is introduced by 
the random signs of the interactions Jij . 

Figure El shows the FORC diagram of the EASG. One 
can see a pronounced ridge along the _ff c -axis together 
with an asymmetric feature at small coercivities. The un- 
derlying details of the EASG FORC diagram have been 
discussed elsewhere Q- 

In Fig. H the SFH of the EASG is shown. Note that 
the asymmetry present in the FORC diagram in Fig. 0] 
is lost. The weight of the asymmetric part of the FORC 
diagram shifts to the ridge at H\, = 0. 

Although some of the features in the FORC diagram 
(Fig. |3J| of the EASG are missing in the corresponding 
SFH (Fig.EJ), the horizontal ridge reminiscent of the un- 
derlying reversal-symmetry of the Hamiltonian || is con- 
served. In particular, by comparing the FORC diagram 
and SFH one can study the effects of frustration on the 
hysteretic behavior of a spin glass. 



CONCLUSIONS 

By re-parameterizing the major hysteresis loop with a 
switching-field histogram we show that for systems with 



FIG. 4: FORC Diagram of the EASG. Note the ridge along 
the Tic-axis. Data for 5000 disorder realizations and N = 50 3 
spins. 



H b 



12 3 



FIG. 5: SFH of the EASG. While the ridge along the tf c -axis 
is qualitatively conserved, the SFH shows drastic differences 
to the FORC diagram presented in Fig. 0] In particular, the 
asymmetry with respect to the horizontal axis is lost. 



no frustration (random-field Ising model) the SFH closely 
resembles the FORC diagram. Small differences can be 
attributed to numerical error in the calculation of an 
FORC diagram, hysteron correlations, or the breakdown 
of the hysteron picture. 

SFHs show more details about the system than the 
major hysteresis loop and are considerably faster to com- 
puter than FORC diagrams. Therefore they are an effi- 
cient alternative in order to study the microscopic distri- 
butions of coercivity and bias of the switching units. 

Because the switching fields of the underlying micro- 
scopic switching units have to be recorded for the com- 
putation of an SFH, the method is in general limited 
to numerical studies of hysteretic systems. Experimental 
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applicability might be possible with synthetic particulate 
samples |l6j where the individual switching units can be 
traced during the magnetic field sweep. 

We also present results on the (frustrated) Edwards- 
Anderson Ising spin glass. We show that there are clear 
differences between the FORC diagram and the SFH be- 
cause SFHs do not take into account multiple switching 
events of the spins, a hallmark of spin glasses. We sug- 
gest these differences can be used to quantify the effects 
of frustration in FORC diagrams. 

The FORC and SFH methods promise to be power- 
ful tools to "fingerprint" hysteretic systems. Still, the 
breadth of information they provides remains to be un- 
derstood fully. 
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